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Abstract.  The  study  demonstrates  that  an  incompati¬ 
bility  between  a  surface  temperature  climatology  and  a 
given  ocean  model,  into  which  the  climatology  is 
assimilated  via  Haney  restoration,  can  cause  model 
ocean  climate  drift  and  interdecadal  oscillations  when 
the  ocean  is  switched  to  a  weaker  restoration.  This  is 
made  using  an  idealized  Atlantic  Ocean  model  driven  by 
thermal  and  wind  forcing  only.  Initially,  the  temperature 
climatology  is  forcefully  assimilated  into  the  model,  and 
an  implied  heat  flux  field  is  diagnosed.  During  this  stage 
any  incompatibility  is  suppressed.  The  restoring  boun¬ 
dary  condition  is  then  switched  to  a  new  forcing 
consisting  of  a  part  of  the  diagnosed  flux  and  a  part  of 
the  restoring  forcing  in  such  a  way  that  at  the  moment 
of  the  switching  the  heat  flux  is  identical  to  that  prior  to 
the  switching.  Under  this  new  forcing  condition,  the 
incompatibility  becomes  manifest,  causing  changes  in 
convection  patterns,  and  producing  drift  and  interde¬ 
cadal  oscillations.  The  mechanisms  are  described. 


1  Introduction 

It  is  an  open  question  as  to  whether  or  not  an  ocean 
general  circulation  model  (OGCM)  with  prescribed 
physics  is  able  to  maintain  a  given  set  of  prescribed 
surface  climatology.  In  areas  characterized  by  sharp 
oceanic  fronts  and  high  convective  activity,  the  OGCM, 
due  to  the  limitation  of  model  dynamics,  cannot 
maintain  the  climatology.  Therefore,  in  these  areas, 
the  OGCM  will  systematically  bias  away  from  the 
climatology  when  it  is  allowed  to  do  so.  We  refer  to  this 
situation  as  an  incompatibility  between  the  climatology 
and  the  OGCM  dynamics.  In  ocean  modeling,  the 
climatology  is  often  assimilated  into  an  OGCM  via  a 
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Haney  (1971)  restoration.  This  procedure  effectively 
prescribes  the  surface  model  solution  to  be  more  or  less 
the  same  as  the  climatology,  disregarding  whether  or  not 
an  incompatibility  exists.  When  such  an  incompatibility 
exists  it  is  suppressed.  Once  it  is  allowed  to  manifest 
itself,  for  example,  switching  from  a  strong  restoration 
to  a  weak  one,  or  to  a  diagnosed  flux  forcing,  what  are 
the  effects  of  the  incompatibility?  This  is  the  central  issue 
of  this  study.  This  is  examined  in  a  pure  thermal 
circulation  system.  We  demonstrate  that  it  will  cause  a 
drift  from  the  state  obtained  by  the  strong  restoration 
and  internal  variabilities. 

Recently,  there  has  been  a  surge  of  research  interest 
in  the  internal  variability  of  the  ocean  thermohaline 
circulation.  Kushnir  (1994)  presented  evidence  showing 
oceanic  variabilities  in  the  North  Atlantic  with  a  30-40 
year  time  scale.  In  a  similar  study,  Deser  and  Blackmon 
(1993)  proposed  that  the  general  warming  trend  in  the 
North  Atlantic  sea  surface  temperature  (SST)  during  the 
1920s-30s  and  the  cooling  in  the  1960s  is  associated  with 
changes  in  the  Gulf  Stream  system.  These  new  results 
emerge  at  a  time  when  results  from  the  Geophysical 
Fluid  Dynamics  Laboratory  (GFDL)  coupled  ocean- 
atmosphere  model  show  the  existence  of  similar  oscil¬ 
lations  in  the  thermohaline  circulation  in  the  North 
Atlantic  with  a  comparable  time  scale  of  40  to  50  years 
(Delworth  et  al. ,  1993).  Deser  and  Blackmon  (1993)  also 
showed  a  variability  of  about  10  years.  Decadal  and 
interdecadal  oscillations  driven  by  a  constant  flux 
forcing  have  been  produced  in  many  models  (e.g., 
Huang  and  Chou,  1994;  Huang,  1994;  Cai  et  al .,  1995; 
Cai  1995).  In  the  Cai  et  al.  (1995)  study  an  idealized 
Atlantic  ocean  is  spun  up  by  restoring  the  model’s 
uppermost  level  temperature  to  a  specified  temperature 
climatology  alone.  When  the  model  reaches  a  steady 
state,  the  forcing  is  switched  to  a  surface  heat  flux 
diagnosed  from  the  steady  state.  This  switch  does  not 
induce  interdecadal  oscillations.  A  small  zonal-redistri¬ 
bution  of  the  diagnosed  heat  flux  then  leads  to  strong 
interdecadal  oscillations  of  the  thermal  circulation.  By 
zonal-redistribution,  it  is  meant  that  the  east-west 
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variation  of  the  surface  heat  flux  is  modified  by  forming 
a  new  flux  field  that  is  a  linear  combination  of  the 
diagnosed  flux  and  its  zonal  mean.  In  this  way,  there  is 
no  addition  or  loss  of  heat  associated  with  the  redistri¬ 
bution,  and  the  oceanic  heat  transport  implied  by  the 
redistributed  flux  field  is  the  same  as  that  implied  by  the 
diagnosed  field.  The  model  SST  anomalies  are  similar  to 
those  observed  in  the  North  Atlantic  (Deser  and 
Blackmon,  1993;  Kushnir,  1994),  and  to  those  of  the 
interdecadal  oscillations  in  the  GFDL  coupled  atmo¬ 
sphere  ocean  model  (Delworth  et  al. ,  1993). 

The  present  study  shows  that  no  heat  redistribution  is 
required  for  the  generation  of  oscillations.  Diagnosed 
heat  fluxes  are  not  necessarily  compatible  with  the 
model  dynamics.  Indeed  if  the  climatology  cannot  be 
maintained  by  the  model  dynamics,  neither  will  be  the 
diagnosed  heat  flux.  In  this  context,  let  us  consider  the 
following  situations.  First,  given  a  surface  climatology 
and  an  OGCM,  an  infinite  numbers  of  states  and 
associated  heat  fluxes  may  be  obtained  via  different 
restoration  time  scales.  Second,  given  an  OGCM  and  a 
restoration  time  scale,  different  flux  fields  may  be 
obtained  with  different  sets  of  climatology.  Finally, 
different  fluxes  may  be  obtained  by  changing  model 
parameters.  It  is  just  likely  that  diagnosed  fluxes  are 
incompatible  with  the  model  dynamics.  We  show  that  an 
incompatibility  is  readily  present  if  the  model  domain  is 
extended  polewards  slightly  from  that  of  Cai  et  al. 
(1995).  This  incompatibility  has  a  significant  impact 
even  with  a  small  change  in  the  thermal  damping 
strength. 

The  remainder  of  this  study  is  divided  into  three 
sections.  In  Sect.  2,  the  model  and  the  thermal  forcing 
are  described.  In  Sect.  3,  model  results  are  presented. 
The  conclusions  and  the  implications  for  coupled 
models  are  given  in  Sect.  4. 


2  The  model  and  the  thermal  forcing  condition 

This  study  employs  the  Pacanowski  et  aV s  (1991) 
version  of  the  Bryan-Cox-Semtner  OGCM,  which  is 
based  on  the  work  of  Bryan  (1969).  The  horizontal  grid 
spacing  is  4°  latitude  by  4°  longitude.  The  model  has  12 
levels  in  the  vertical,  at  depths  listed  in  Table  1.  Values 
assigned  to  the  various  model  parameters  are  listed  in 
Table  2.  These  are  the  same  as  Cai  et  al.  (1995),  Cai  and 
Godfrey  (1995),  or  Cai  (1995).  The  model  domain  is  a 
flat-bottomed,  two-hemisphere  Atlantic  basin.  It  has  a 
width  of  60°  and  a  length  of  144°  extending  from  72°  S  to 
72°N.  This  extends  from  the  model  of  Cai  et  al.  (1995)  in 
the  north-south  direction  by  one  model  grid  in  each 
model  hemisphere.  The  Southern  Hemisphere  includes  a 
modeled  Antarctic  Circumpolar  Current  (ACC)  passage 
from  44° S  to  60° S,  and  a  sill  2350  m  deep  in  the  model’s 
Drake  Passage.  This  sill  is  intended  to  provide  a  realistic 
ACC  transport  through  the  setup  of  a  bottom  pressure 
difference  between  the  two  sides  of  the  passage  (Hol¬ 
land,  1973;  Gill  and  Bryan,  1971;  Cai,  1994).  Otherwise 
the  flat  bottom  would  lead  to  an  unrealistically  large 
ACC,  and  an  ACC  has  to  be  specified  (Cai  et  al .,  1995). 


Table  1.  Distribution  of  vertical  levels 


Level 

Thickness  (m) 

Depth  of  (T,S)  (m) 

1 

25.0 

12.5 

2 

25.0 

37.5 

3 

40.0 

70.0 

4 

70.0 

125.0 

5 

110.0 

215.0 

6 

200.0 

370.0 

7 

330.0 

635.0 

8 

450.0 

1025.0 

9 

650.0 

1575.0 

10 

900.0 

2350.0 

11 

900.0 

3250.0 

12 

900.0 

4150.0 

Table  2.  Values  of  model  coefficients 


Parameter 

Symbol 

Value 

Horizontal  diffusivity 

ATii 

1  x  103  m2  s-1 

Horizontal  viscosity 

^4mh 

3  x  105  m2  s-1 

Vertical  diffusivity 

Ajy 

1  x  10-4  m2  s-1 

Vertical  viscosity 

A  MV 

1  x  10-4  m2  s-1 

Another  small  modification  from  that  of  Cai  et  al. 
(1995)  is  the  inclusion  of  a  model  Weddell  Sea  to  allow 
an  Antarctic  Bottom  Water  Formation  cell. 

A  simple  way  in  which  the  incompatibility  manifests 
is  through  a  change  in  thermal  damping.  During  an 
ocean  spinup,  in  order  to  best  represent  the  current 
ocean  climate,  strong  restorations  are  usually  used. 
Upon  switching  to  an  active  atmosphere,  the  thermal 
damping  weakens,  therefore  SST  in  a  coupled  environ¬ 
ment  is  more  free  to  evolve,  and  any  incompatibility  is 
manifest. 

In  the  present  study  we  use  a  simple  thermal 
boundary  condition.  This  is  diagnosed  from  the  follow¬ 
ing  procedure.  Suppose  a  restoring  spinup  is  carried  out 
in  which  the  model  SST  is  restored  to  Ta ,  and  the  steady 
state  solution  has  an  SST  field  Ts  and  a  heat  flux  field 
Fspin.  They  satisfy 

Fspin=KH(Ta-Ts).  (1) 

Here  KH  is  the  thermal  damping  rate  during  the  spinup. 
The  heat  flux  Qao  is  switched  to 

Qao  =  (1  -  a )Fspm  +  V-KH(Ta  -  71),  (a  <  1).  (2) 

Here  a KH  is  the  new  damping  rate  and  T\  is  the  time- 
dependent  temperature  of  the  uppermost  level.  At  the 
moment  of  switching  forcing,  the  heat  flux  is  exactly  the 
same  as  that  implied  by  the  steady-state  of  the  spinup. 

The  physics  behind  this  treatment  of  surface  thermal 
forcing  can  be  realized  from  a  simple  coupled  system 
consisting  of  a  simple  atmosphere  and  a  mixed  layer 
ocean  (Cai  et  al .,  1995;  Cai  and  Chu,  1996).  The  heat 
balance  for  the  atmosphere  is 

0  =  -KH(Ta  -  71)  -  K'rTa  +  Qa,  (3) 

and  for  the  mixed  layer  ocean  is 
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C°~dt^  KH(Ta**  T\)  +  Q0.  (4) 

Equation  (3)  expresses  a  balance  between  heat  transfer 
from  the  ocean  to  the  atmosphere  given  by  the 
—KH(Ta  —  T\)  term,  heat  loss  to  space  given  by  — K'rTa , 
and  heat  sources  due  to  solar  input  and  the  divergence 
of  the  atmospheric  heat  transport  given  by  Qa.  Here,  K'r 
is  the  atmospheric  radiative  feedback  parameter,  CQ  is 
the  specific  heat  capacity  of  the  mixed  layer  ocean,  and 
Q0  is  the  oceanic  heat  source. 

Defining  a  temperature 

Tr  =  QJK,  (5) 

we  can  write  Eq.  (3)  as: 

0  =  -Kff(Ta  -  Ti)  +  K'r(Tr  -  Ta).  (6) 

From  Eq.  (6),  expressing  Ta  in  terms  of  other  quantities 
and  substituting  into  Eq.  (4),  we  have: 

C0^  =  Kr(Tr-Tl)  +  Q0,  (7) 

where 


Kr  = 


KlrK'r 

(Kh  +  K>r) ' 


(8) 


This  means  that  when  the  mixed  layer  ocean  is  coupled 
to  the  simple  atmosphere  described  by  Eq.  (3),  its  heat 
balance  can  be  described  by  Eq.  (7).  Comparing  Eqs.  (4) 
and  (7),  we  see  that  the  heat  exchange  term 


KH(Ta  —  T\)  =  Kr(Tr  —  T\).  (9) 

Defining  a  ratio 

a  =  Kr/KH  (a<l),  (10) 


which  measures  the  change  in  thermal  damping,  then  we 
can  rewrite  the  heat  exchange  term  of  the  simple 
coupling,  using  Eqs.  (1)  and  (9),  as 


Kr(Tr  -  Tx)  =Kr(Tr  -  Ts)-Kr(Ta  -  Ts)  +  Kr(Ta  -  T\) 

=  (1  -  <x)Fspin  +  uKH(Ta  -  7i).  (11) 

Thus,  the  new  forcing,  which  consists  of  a  part  of  the 
constant  diagnosed  heat  flux  plus  a  weaker  restoring 
forcing,  is  equivalent  to  the  simple  coupling.  Note  that 
when  switching  from  a  spinup  using  Haney  restoration 
as  a  thermal  boundary  condition  to  one  provided  with  a 
fully  coupled  model,  a  changes  from  1  to  a  value  ranging 


from  0.1  to  0.2.  In  the  present  study,  we  show  that  an 
incompatibility  can  manifest  even  with  a  small  change 
in  a. 


3  Model  results 

3.1  Model  spinup 

Initially,  the  OGCM  is  spun  up  by  restoring  the 
uppermost  level  temperature  a  zonally  uniform  temper¬ 
ature  given  by  the  profile  shown  in  Fig.  la,  which 
approximates  the  observed  zonally  averaged  field.  This 
run  is  referred  to  as  R15  (see  Table  3).  The  profile  has 


Fig.  la.  Surface  relaxation  temperature  (in  °C),  and  b  overturning 
stream  function  (in  Sverdrup)  from  R15 


Table  3.  Details  of  model  ex¬ 
periments  described  in  Sect.  3. 
The  second  column  indicates 

Run 

Initial 

condition 

Model 

domain 

Forcing 

Persistent 

oscillations 

initial  condition.  The  fourth 

R15 

Rest 

72°S  to  72°N 

Restoring  (15  days) 

None 

column  shows  the  forcing  con- 

C30 

R15 

72°S  to  72°N 

New  forcing  (30  days) 

None 

dition  and  the  damping  time 

C45 

R15 

72°S  to  72°N 

New  forcing  (45  days) 

27.1  y 

scale  is  indicated  in  the  brack¬ 

C45r 

R15 

72°S  to  72°N 

Restoring  (15  days)  68°N 

None 

ets.  The  fifth  column  indicates 

northward  and 

whether  or  not  persistent  oscil¬ 

New  forcing  (45  days) 

None 

lations  occur.  If  yes,  the  period 

elsewhere 

is  given  in  years  (y) 

R15CGZ 

Rest 

68°S  to  68°N 

Restoring  (15  days) 

None 

CGZ300 

R15CGZ 

68°S  to  68°N 

New  forcing  (300  days) 

None 

R45 

R15 

72°S  to  72°N 

Restoring  (45  days) 

24.7  y 
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the  same  analytical  expression  as  in  Cai  et  al.  (1995)  and 
the  expression  applies  to  the  extended  model  domain, 
that  is,  from  68°  to  72°.  The  e-folding  restoring  time  is 
15  days.  This  corresponds  to  a  thermal  damping  KH  of 
79Wm-2oC-1.  The  model  is  also  subject  to  the  wind 
stress  used  by  Bryan  (1987).  Only  during  the  spinup  are 
the  lower  level  acceleration  techniques  of  Bryan  (1984) 
used.  After  about  4000  surface  years  (32  000  years  at 
bottom  level)  of  integration,  the  spinup  reaches  a 
statistically  steady-state.  The  overturning  circulation  is 
shown  in  Fig.  lb,  and  is  qualitatively  similar  to  that  in 
the  real  Atlantic.  It  features  a  northern  sinking  cell;  a 
southern  sinking  cell;  and  an  Antarctic  Circumpolar 
Current  Convective  cell.  The  model  ACC  reaches 
104  Sv. 

The  convection  pattern  is  shown  in  Fig.  2.  During  the 
spinup,  SST  is  forcefully  “nudged”  to  the  prescribed 
climatology,  regardless  of  whether  or  not  the  internal 
dynamics  can  maintain  it.  This  nudging  in  effect 
prescribes  the  convection  pattern,  and  the  heat  flux  for 
its  maintenance.  Any  incompatibility  between  the 
assimilated  climatology  and  the  model  dynamics  is 
suppressed. 


3.2  Drift  and  variability  under  the  new  forcing  condition 

In  the  following  two  experiments,  the  thermal  damping 
is  switched  to  the  new  forcing  described  by  Eq.  (2),  with 
a  damping  rate  equivalent  to  30  (a  =  0.5)  and  45  days 
(a  =  0.333),  respectively.  These  experiments  are  referred 
to  as  C30  and  C45  (see  Table  3).  The  change  in  thermal 
damping  is  small,  but  this  is  intended  to  show  the  high 


1  23456789  10 

Convection  pattern,  R15 


0  30 E  60 E 

Fig.  2.  Convection  penetration  depth  for  the  run  R15.  The  penetra¬ 
tion  depth  is  shown  in  terms  of  model  levels 


sensitivity  to  this  parameter  in  the  presence  of  an 
incompatibility. 

If  the  model  dynamics  and  the  assimilated  climatol¬ 
ogy  are  fully  compatible  with  each  other,  then  under  the 
new  forcing  the  ocean  state  should  not  change,  because 
the  new  forcing  provides  an  identical  air-sea  heat  flux  at 
the  moment  of  the  switching.  In  the  presence  of  an 
incompatibility,  however,  the  ocean  undergoes  signifi¬ 
cant  changes. 

Figure  3  shows  the  time  series  of  the  overturning 
sampled  at  a  location  of  (60°N,  2350  m).  Unless 
otherwise  stated,  single-point  time  series  of  overturning 
are  all  taken  at  this  location.  Note  that  this  does  not 
coincide  with  the  location  where  the  overturning  under 
the  new  forcing  is  a  maximum,  nor  with  the  location 
where  the  amplitude  of  the  oscillation  is  a  maximum 


Fig.  3a, b.  Time  evolution  of  the  northern  overturning  (in  Sverdrup) 
at  a  location  of  (64°  N,  2350  m)  for  a  C30  and  b  C45 
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under  the  flux  corrected  condition.  Several  features 
emerge.  First,  in  run  C30  (Fig.  3a),  the  steady-state 
drifts  from  that  of  R15  with  some  initial  oscillations. 
The  time  scale  for  the  drift  is  several  thousands  of  years. 
Second,  in  run  C45  (Fig.  3b),  besides  a  similar  drift 
regular  thermal  oscillations  with  a  period  of  27. 1  y  persist. 

To  understand  the  cause  for  these  features,  we  focus 
on  the  results  of  C45,  since  these  features  are  more 
profound  in  this  case.  Figure  4  shows  the  convection 
patterns  at  a  quarter  period  of  the  steady  oscillation.  We 
see  that  as  the  overturning  oscillates,  so  does  the 


Fig.  4a-d.  The  same  as  in  Fig.  2  but  for  run  C45  at  a  quarter  period 
interval.  Panels  b  and  d  correspond  to  times  when  the  overturning  at 
the  location  of  (64°  N,  2350  m)  is  at  a  minimum  and  maximum 
respectively 


convection  pattern,  especially  in  the  regions  60°N 
northward  and  off  the  northeastern  boundary.  These 
changes  in  the  convection  patterns  produce  significant 
drifts  of  temperatures,  especially  SSTs,  which  are  shown 
in  Fig.  5,  again  at  a  quarter  period.  We  see  that  drifts  at 
peak  and  trough  are  not  antisymmetric.  It  suggests  that 
the  mean  state  has  also  drifted  from  the  spinup  state. 
Figure  6a, b  shows  the  drift  of  SST  and  the  temperature 
at  2350  m  of  the  mean  state  averaged  over  10  steady 
cycles.  At  the  surface  the  largest  drift  occurs  in  the 
convection  regions,  and  at  depth  the  ocean  is  generally 
cooled.  An  examination  reveals  that  the  cooling  is 
associated  with  a  decrease  in  the  northern  high-latitude 
sinking  cell. 


a  SST  drift  at  mid.  peak 
to  trough 


60N 


30N 


0 


30S 


60S 


0  30E  60E 

c  SST  drift  at  mid.  trough 
to  peak 


0  30E  60E 


b  SST  drift  at  trough 


Fig.  5a-d.  The  same  as  Fig.  4  but  SST  drift  (in  °C)  from  the  final  state 
of  R15 
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a  Mean  SST  drift  b  Mean  temp,  drift  at  level  10 
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Fig.  6a, b.  Drifts  of  the  mean  state  of  C45  from  the  steady-state  of 
R15,  a  SST  and  b  temperature  at  2350  m.  The  mean  state  is  obtained 
by  averaging  over  10  steady  cycles 

These  features  are  associated  with  instability  of 
convections  of  the  spinup  state.  In  the  presence  of  an 
incompatibility,  the  equilibrium  solution  of  the  spinup 
via  a  strong  relaxation  to  a  given  climatology  is  not  a 
perfect  steady-state  solution  of  the  model  equations.  In 
areas  characterized  by  sharp  oceanic  fronts  and  high 
convective  activity,  the  model,  due  to  the  limitation  of 
model  dynamics,  cannot  maintain  the  climatology. 
Therefore,  in  these  areas,  the  model  is  trying  to 
systematically  move  away  from  the  climatology.  How¬ 
ever  under  a  strong  restoration,  the  surface  solution  is 
strongly  constrained.  Consequently,  globally  averaged 
implied  surface  heat  and  salinity  fluxes  (figures  not 
shown)  tend  to  oscillate  about  some  mean  value,  and 
convections  are  intermittent.  When  the  restoration 
weakens,  the  surface  model  field  is  allowed  to  evolve 
so  that  it  is  more  compatible  with  the  dynamics.  The 
pattern  of  convections  change,  and  a  dynamical  adjust¬ 
ment  ensues.  The  features  of  drift  and  oscillation  are  a 
manifestation  of  the  adjustment  process  suppressed 
during  the  spinup.  Since  the  incompatibility  is  present 
permanently,  as  the  dynamical  adjustment  proceeds,  the 
system  is  constantly  seeking  a  new  convection  pattern, 
giving  rise  to  the  oscillation. 

The  oscillation  in  the  present  study  is  associated  with 
a  dynamic  adjustment  process  in  a  form  of  waves  that 
propagate  along  the  model  basin.  Detailed  characteris¬ 
tics  are  discussed  in  Cai  and  Chu  (1997).  In  short,  the 
oscillation  is  associated  with  generations  and  propaga¬ 
tions  of  long  Planetary  waves,  that  is,  westward  propo- 
gating  Rossby  waves.  This  is  in  sharp  contrast  to  results 
of  Winton  (1995)  and  Greatbatch  and  Peterson  (1996) 
studies.  They  show  that  oscillations  in  their  models  are 
associated  with  propagations  of  waves.  But  they  believe 
that  the  waves  in  their  models  are  Kelvin  waves, 


although  strictly  speaking  Kelvin  waves  are  not  resolved 
in  OGCMs  of  their  resolutions. 


3.3  Comparison  with  the  results  of  Cai  et  al.  ( 1995) 

In  order  to  locate  where  the  incompatibility  arises,  we 
compare  the  present  results  with  those  of  the  Cai  et  al. 
(1995)  study.  Although  a  similar  drift  is  produced  in  Cai 
et  al.  (1995),  no  oscillation  is  induced  when  the  forcing  is 
switched  to  one  with  a  similar  damping  rate.  As 
commented  upon  earlier,  the  model  here  extends  4° 
poleward  (in  both  hemispheres).  The  fact  that  the 
convection  pattern  oscillates  mainly  in  this  northern 
extended  region  suggests  that  most  of  the  incompatibil¬ 
ity  arises  from  this  region.  To  test  this  point,  two  more 
experiments  are  carried  out.  Run  CGZ300  is  conducted 
in  a  new  model  with  a  north-south  extension  corre¬ 
sponding  exactly  to  that  in  Cai  et  al.  (1995),  that  is,  from 
68° S  to  68°N,  and  under  a  new  damping  corresponding 
to  300  days.  The  model  experimental  procedure  and 
model  coefficients  are  the  same  as  described  already 
involving  an  initial  spinup  from  rest  (R15CGZ,  see 
Table  3).  In  this  experiment,  switching  to  a  new 
damping  produces  no  oscillation,  even  with  a  large 
change  in  a,  i.e.,  from  1  to  0.05.  Figure  7a  shows  the 
evolution  of  the  overturning.  Although  some  weak 
initial  oscillations  are  evident,  they  die  out  soon.  Run 
C45r  is  identical  to  C45,  except  that  in  the  northern 
extended  region  (that  is,  the  row  of  northernmost  grids) 
the  restoring  forcing  (hence  the  ‘r’)  is  the  same  as  in  the 
spinup.  Figure  7b  shows  the  time  series  of  the  over¬ 
turning.  Again  no  oscillation  is  produced.  These  two 
experiments  indicate  that  the  incompatibility  responsible 
for  the  oscillation  arises  from  the  extended  region.  These 
also  suggest  that  the  oscillatory  behavior  of  the  convec¬ 
tion  pattern  in  other  northern  high-latitudes  in  C45  is 
induced  by  the  incompatibility  in  the  extended  region. 

3.4  Results  under  a  weak  restoration 
with  no  added  flux  term 

We  have  argued  that  the  incompatibility  between  the 
climatology  and  the  internal  model  dynamics  is  respon¬ 
sible  for  the  oscillation  in  C45.  Such  a  notion  implies  that 
if  the  incompatibility  is  not  suppressed  by  the  strong 
restoration  as  in  R15,  it  will  manifest  itself  and  generate 
similar  oscillations.  To  this  end,  R45  is  carried  out,  in 
which  the  temperature  of  uppermost  level  is  relaxed  to 
the  climatology  (Fig.  la)  with  a  restoring  time  of  45  days 
(see  Table  3).  This  experiment  starts  from  the  steady-state 
of  R15,  therefore,  the  only  difference  from  C45  is  that  no 
added  flux  is  employed.  At  the  moment  of  changing  the 
thermal  damping  rate,  the  surface  heat  flux  is  suddenly 
and  uniformly  reduced  to  one  third  of  that  prior  to  the 
change.  The  reduction  of  the  heat  flux  initially  shocks  the 
system  somewhat,  but  upon  recovery  the  incompatibility- 
induced  oscillation  appears.  The  period  is  24.7  (Fig.  8). 

Under  this  weak  restoration,  the  modeled  SSTs  drift 
away  from  the  steady-state  of  R15o.  Figure  9  shows  the 
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Fig.  7a, b.  Time  evolution  of  the  northern  overturning  (in  Sverdrup) 
at  a  location  of  (64°  N,  23507  m)  for  a  CGZ300  and  b  C45r.  See  text 
for  details  about  these  runs 
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Fig.  8.  Time  evolution  of  the  northern  overturning  (in  Sverdrup)  at  a 
location  of  (64°N,  2350  m)  for  R45 
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drift  of  the  mean  state  averaged  over  10  steady  cycles  of 
the  temperature  at  surface  level  and  at  level  10.  The 
surface  drift  (Fig.  9a)  is  much  larger  in  comparison  with 
that  of  C45.  In  particular,  a  substantial  cooling  at  tropics 
and  a  warming  at  high  latitudes  take  place,  in  contrast  to 
the  drift  under  an  added  flux.  The  warmer  SST  is  then 
transmitted  to  the  deep  ocean  by  convection,  leading  to  a 
general  rise  in  temperature  at  depth  (Fig. 9b). 

4  Discussion  and  summary 

An  OGCM  is  usually  spun  up  via  Haney  restoration  to 
the  present  day  climatology  to  an  “equilibrium”  state. 
The  climatology  is  forcefully  nudged  to  the  model  and 
any  incompatibility  between  the  climatology  and  the 


dynamics  is  suppressed.  In  this  study  we  demonstrate 
that  when  under  a  weaker  thermal  damping  environ¬ 
ment  the  incompatibility  manifests,  generating  drift 
(from  the  spinup)  and  internal  variability.  The  process 
is  associated  with  instability  of  convection,  which  occurs 
when  the  surface  forcing  changes  from  a  strong 
restoration  to  a  weaker  one.  When  the  thermal  damping 
is  weak  enough,  oscillations  occur  as  a  result  of  a 
dynamic  adjustment  process  to  the  permanent  presence 
of  the  incompatibility.  The  incompatibility  constantly 
generates  perturbations,  which  propagate  as  boundary 
waves  along  the  weakly  stratified  northern  boundary. 

It  is  striking  that  the  generation  of  oscillations  is 
extremely  sensitive  to  the  incompatibility  at  northern¬ 
most  grids.  A  question  arises  as  to  how  sensitive  this 
result  is  to  the  surface  climatology  used  to  drive  the 
model.  To  answer  this  question,  we  have  carried  out 
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a  Mean  SST  drift  b  Mean  temp,  drift  at  level  10 
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Fig.  9a, b.  Drifts  of  the  mean  state  of  R45  from  the  steady  state  of 
R15,  a  SST  and  b  temperature  at  2350  m.  The  mean  state  is  obtained 
by  averaging  over  10  steady  cycles 

numerous  experiments  varying  the  north-south  gradient 
of  the  restoring  temperature.  In  all  experiments  drift  and 
oscillation  occur  when  the  restoring  time  scale  is 
changed  from  15  days  to  45  days,  suggesting  that  it 
may  be  difficult  to  achieve  a  compatibility.  Given  the 
fact  that  instability  of  a  convection  pattern  plays  a 
significant  role  in  the  generation  of  drift  and  oscillation, 
it  is  highly  relevant  to  ask  how  sensitive  the  model  result 
is  to  the  convection  scheme  used.  To  address  this 
question,  we  have  repeated  run  R15  and  C45  using  the 
convection  scheme  of  Rahmstorf  (1993).  The  advantage 
of  the  scheme  over  Cox’s  (1987)  is  that  it  achieves  a 
complete  mixing  when  vertical  instability  occurs.  We 
found  that  upon  switching  to  a  weaker  thermal  dam¬ 
ping,  the  model  behaves  in  a  similar  way  to  that  seen  in 
the  original  R15o  and  R45. 

That  oscillations  may  be  driven  by  incompatibility 
between  surface  forcing  and  ocean  dynamics  may  have 
significant  implications.  On  the  one  hand,  in  reality, 
there  is  no  dynamical  reason  why  forcing  fields  should 
always  be  compatible  with  the  oceanic  internal  dynam¬ 
ics.  For  example,  an  atmospheric  heat  flux  is  influenced 
by  many  atmospheric  processes  which  are  weakly 
dependent  on  or  even  completely  independent  from 
the  ocean  (Cai  et  al. ,  1995;  Cai,  1995).  Thus  there  may 
exist  some  incompatibility.  In  this  sense,  the  incompat¬ 
ibility  induced  modeled  oscillation  may  simulate  a 
realistic  process.  On  the  other  hand,  there  are  many 
factors  that  can  cause  an  incompatibility  during  a 
spinup.  These  include  inaccuracy  in  the  surface  clima¬ 
tology,  deficiency  in  model  dynamic  processes,  either 
absent,  poorly  parametrized,  or  not  accommodated  by 
the  model  resolution.  For  example,  at  polar  latitudes, 
sea-ice  dynamics,  which  is  not  included  in  the  present 
model,  may  be  essential.  Another  factor  that  may  cause 


incompatibility  is  that  under  polar  conditions,  especially 
in  the  presence  of  sea-ice,  observations  are  often 
difficult,  and  hence  a  climatology  for  the  region  may 
be  poorly  constructed. 

Pre-coupling  spinup  of  an  OGCM  via  Haney  resto¬ 
ration  to  the  present  day  climatology  is  commonly 
carried  out  to  provide  initial  conditions,  and  to  con¬ 
struct  flux  adjustment  fields  (Sausen  et  al .,  1988)  for 
coupled  models  (e.g.,  Manabe  and  Stouffer,  1993,  1994; 
Gordon  and  O’Farrell,  1997).  Nakamura  et  al.  (1994) 
show  that  flux  adjustments  do  not  correct  errors  in 
climate  sensitivity.  Marotzke  and  Stone  (1995)  confirm 
this  conclusion.  Schneider  (1996)  points  out  that  the 
manner  in  which  flux  adjustment  is  currently  applied  is 
inappropriate  in  that  the  adjustment  does  not  eliminate 
model  errors.  Our  study  indicates  that  in  the  presence  of 
an  incompatibility,  a  flux  adjustment,  although  designed 
to  remove  the  systematic  mismatch  between  the  atmo¬ 
sphere  and  the  ocean  components,  has  “side  effects.” 
One  of  them  is  that  it  will  carry  the  incompatibility  to  a 
coupled  model,  and  affect  the  coupled  model  behavior. 
It  seems  that  only  when  incompatibility  is  ultimately 
eliminated  from  the  spinup  stage,  can  flux  adjustments 
truly  perform  their  function  without  this  side  effect.  To 
achieve  this  end,  a  good  climatology  and  parametriza- 
tion  of  model  physics  are  essential. 
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